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Abstract 

We solve the Multichannel Kondo model with channel anisotropy using the 

Bethe Ansatz method. The model generates energy scales, characterizing the 

neighborhoods of the various infrared fixed points, reflecting the structure of 

the symmetry breaking in the channel sector. The nature of these fixed points 

also depends on the magnitude of the impurity spin S. We present a detailed 

discussion for the two channel case and point out some new non-Fermi liquid 

behavior. 
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Non-Fermi liquid (NFL) behavior observed in some Ce and U alloys [|l| has stimulated 
an intense study of many impurity models; in particular, of the multichannel Kondo model 

/ 

= i/o + 2 E E Jm^i,™(0)o-,,^f,,™(0) ■ S, (1) 

a,b m=l 

where Hq = —iJ2a;mI dxip\^{x)dx'>pa,m{x). Here the field ilJa,m describes electrons with spin 
index a = ±1/2 and orbital channel [flavor) index m, m = 1.../, and we chose to set vp = 1. 
The operator S represents the spin-S* impurity localized at x = 0. 

The infrared behavior of the model depends on the structure of the flavor sector [Q. In 
particular, a new behavior appears when / > 25", with the overscreened system exhibiting 
NFL physics. 

An exact solution was found for the isotropic case, Jm = J • The solution yields the 
spectrum and the thermodynamics for any temperature and magnetic field. In particular, the 
residual entropy and the critical exponents governing the low-energy physics were calculated. 
The thermodynamics was studied numerically in 

The neighborhood of the isotropic fixed point was further studied by means of Conformal 
Field Theory g| and Bosonization methods [0, yielding the long distance asymptotics of 
the correlation functions with critical exponents that are same as those characterizing the 
thermodynamic functions. 

In this paper we present the solution of the channel-anisotropic model. We shall discuss 
in detail the case of two channels and arbitrary spin S, and briefly outline the generalization 
for more flavors. 

We shall find that in the two channel case the model generates two scales, Tj and T^, which 
we shall later interpret as associated with two fixed points, the isotropic and anisotropic, 
respectively. The scales are explicitly given (in our cutoff scheme) by = De •-'i , Ta = 
D cos is the number of electrons in each channel, L is 

the length of the system, and Ji < J2. In the scaling limit D is taken to infinity with the 
(bare) couplings' dependence on D chosen so as to keep the scales finite. The functional 
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dependence of the scales on the couphng constants Ji, J2, is not universal and may change 
with the cutoff procedure. However, the dependence of physical quantities on the scales is 
universal. The ratio A = T^/Tj is the physical measure of the anisotropy. 

The presence of flavor allows up to / electrons to interact simultaneously with the impu- 
rity. Therefore, spin composites of electrons form irrespectively of the degree of anisotropy, 
and their binding energy is set by the smallest of the couplings. The hamiltonian must be 
regularized with care to allow the formation of the composites while maintaining integra- 
bility. We choose the regularization scheme used in f^, leading to the (regularized) first 
quantized form of the hamiltonian (|l]), 

2N 

h = E {-^d, - A-\d,f + 2J,5{x,)a, ■ S) . (2) 
i=i 

Here Jj is either Ji or J2- The limit A ^ 00 will be taken only after determining the 
eigenvalues. This limiting procedure, the fusion |^, leads to the formation of the spin 
composites. We shall find that the binding energy will be well above the spin energy scales, 
and the flavor excitations will disappear from the low-energy spectrum. 

The eigenfunctions of are combinations of plane waves with pseudomomenta {kj,j = 
1, . . . , 2A^} and amplitudes ^^^;;;^^q, depending on the electron spin and flavor indices 
aj,mj, and the impurity spin index a. The energy eigenvalues in terms of the pseudomo- 
menta are E = I]j=i + %/A), while the amplitudes are determined from the two body 
S-matrices, to which we now turn. 

The impurity-electron S-matrix is derived from (|^), 

_ A, + 1-.J,(^-S + |) 
^jo , 1-5 j 

Aj + 1 — tJj 

where Aj = kj/A. Since the interactions are flavor preserving, the S-matrix (j^) has only a 
non-trivial term in spin space. 

All electrons are right-movers and as there is no direct interaction between them the elec- 
tronic states are infinitely degenerate away from the impurity. This permits the introduction 
of an arbitrary electron-electron S-matrix Sji into the definition of the wave-functions. To 
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construct a basis that will manifest the integrability of the model the matrices Sji must be 
chosen so as to satisfy the Yang-Baxter factorization equations, 

SjoSioSji = SjiSioSjo, SjkSikSji = SjiSikSjk- (4) 

The solution for Sji can be given as a direct product of spin and flavor terms, Sji = S'jf*" ® 
gfiavor ^ with each term satisfying (H) separately. The spin component of Sji is given by 

where Xj = Xj/ Jj and Pji = \ {cTj ■ <Ji + 1). The form of @ reflects the SU{2)spin invariance 
of the model. The flavor component of Sji reflects the breaking of the SU (2) flavor symmetry 
to a residual U{1) when the anisotropy is present, and is given by 



^flavor ^ 1 ^ n-'^'^-N , • ^ ®r, + ry® Ty} 

^ smh(K(Aj — A/) + lu) 

1 1 1 sinh(fi;(Ai — A^)) 

+ -{1 + ® rj + \^ ^ ^ — — {1 -T^(g) rj. 

2 2 sinh(K(Aj — A^) + iu) 

Here, {r} are the Pauli matrices, and k and z/ are functions of the couplings. Denoting 
fi = u/hi, we shall see that /i is related to the binding energy of the composites. 

Given the set of consistent S-matrices we may derive the Bethe Ansatz equations deter- 
mining the allowed pseudomomenta kj and hence the spectrum. Introducing the auxiliary 
variables {uj-y} for the flavor sector, and {x-y} ioi the spin sector, we find: 



M 



Xs-X,-j- + ^,M, sinh(z/(u;, - + |)) 



e 



Yi - - '-f i Yi 

5=1 Xs - Xj - J- - ^ sinh(z/(to'5 - ^ - f )) 
sinh(z/K - ^ + t)) Ml smhiuju, - cos + |)) 
j=i sinh(z/(cu^ - 7f - f )) -5=1 sinh(i/(w^ - a;^ - f)) ' 
^ X,-Xs + ^ _ X, + ^S '^ X,-X,-j- + ^, 

5=1 ~ XS — i X-y ~ '^S -^^ X-y ~ Xj — J- ~ 2 



describing the full content of the model. 

The ground state and low lying energy excitations reside in the flavor-singlet sector where 
the solutions are in the form of double 2-strings, 



/i , I 



Jl 



(6) 



The structure of the double string solution reflects the flavor symmetry breaking. In the 
isotropic limit the two strings coalesce leading to the SU (2) flavor degeneracy and the NFL- 
behavior. 

The energy associated with the double 2-string is £,5 = 2^u!s (1 + /icu^/A) — /i^A/2, leading 
to the 1 -string hypothesis for the {uj^}, namely us = ps/{iiK) |]^. We substitute the double 2- 
string solution (P) in the eigenvalue equations and take the A — >■ (X) limit. As explained, this 
is the fusion process in which the flavor degrees of freedom are removed from the low energy 

JV 

effective spin problem. We flnd that the spin contribution to the energy is E = I]/=i 4p5, 
and the spin degrees of freedom are described by the following system of fused Bethe Ansatz 
equations: 
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-n 



5=1 Xj Xs 



— I 



X7 - 



X 



X 



X7 



J2 



.^7 J2 2(-'^"'" 



^7 



X7 



Jl 



0) 



X7-7r-t(i + ^)y Vx7-7r-|(i-v^) 

with (f) = /i/ J2 and = fi/Ji. As the binding energy is set by the weakest coupling, fi is 
related to Ji, and in the scaling limit we flnd ^ = Ji, = J1/J2, f = 1. Note that the 
equations reduce to the isotropic equations for = J1/J2 = 1 and to the 1-channel Kondo 
equations for = 0, Ji = 0. 

We now solve the equations, identify the ground state and excitations, and summing 
over the latter, derive the free energy. The main results are: 

{i) The solutions of the equations are of the form given by the string hypothesis valid in 
the thermodynamic limit, Xs''^ = + + 1 ~ 2fc), k = 1, . . . ,n, Xs real. 

(ii) The ground state is composed of x 1" ^ind 2-strings, interpolating between the 
isotropic model (where = 1 and the ground state built of 2-strings leading to NFL physics). 



and the single channel model (where = and the ground state is built of 1-strings de- 
scribing a Fermi liquid (FL)). 

(Hi) The impurity free energy for impurity spin-^*, anisotropy A, temperature T and 
magnetic field h is 

T /-oo ^^ln(l + r725(e,_|_)) 

T - 



■J poo J 

F\i^,S-T,h)^-—j_Ji- 



-oo cosh(^ + ln^) 

The function 7725(C) h/T) belongs to the set of functions {rjn} satisfying the following set of 
coupled integral equations, here written in the scaling limit {D — > 00, keeping Tj, fixed), 

In rji = -2Ae^ + G ln(l + 772) 

In 772 ^ -e^ + G ln(l + 771) + G ln(l + 773) 

ln77„ = G'ln(l + 77„_i) + G'ln(l + 77„+i), n>2, 

with boundary conditions lim„^oo([?^+l] ln(l-|-77„) — [ti] ln(l-|-77„+i)) = —2ixh/T. The integral 

operators [a] and G are defined by the kernels a/(27r(((^' — ^)'^ + (f )^)) and 1/(2 cosh(7r(C — 
C'))), respectively. 

We now elaborate on {iii) . We studied analytically the equations in the hmits of small 
and large anisotropy, both for = 1/2 and for S = 1, and deduced the infrared behavior of 
the thermodynamic functions and its dependence on the anisotropy A: 

The S — 1/2 impurity contribution to the free energy at low temperatures is given by 



^ A < 1 



-f In 2 + I (7,. + a;,. (A)^) f In J + 0(f ) + 0(Af ), « T « 

. -Ki - i) h.h^-.h (^)' + ^(^^-^)) S' ^ « ^' ^ » 1 

where 7^ 1 = 7r/6, cu^ 1 = I/tt are the coefficients for the specific heat and susceptibility 
of the S — 1/2 single channel Kondo model, and 7^1 =1/4, a;-_i = 2/7r^ are those of the 
S — 1/2 two-channel isotropic model, and 5 is a constant close to 1. 

The free energy yields the susceptibility and the specific heat. We find, for any A > 0, 
a linear specific heat and a temperature independent susceptibility. Therefore, any small 
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amount of anisotropy moves the system away from the isotropic fixed point. In particular, 
the I In 2 contribution to the T = entropy that appears in the isotropic case is no longer 
present. However, for small anisotropy the behavior of the system in the temperature range 
Ta = ATi < T < Tj is identical to that of the isotropic case, and the isotropic NFL behavior 
reemerges upon setting A = 0. For large anisotropy, on the other hand, the leading terms of 
the thermodynamic functions are those of the S = 1/2 single channel model with Ta playing 
the role of the Kondo temperature. 
Computing the ratio R, 



Ri 

2 



3' 



A = 

-8AlnA, 0<A<1 
2 + C(Ae-^), A > 1 



we again conclude that turning on the anisotropy destroys the NFL behavior of the isotropic 
model, while the A = cxd single channel FL is only weakly modified upon reducing the 
anisotropy. 

We now turn to = 1. The impurity free energy and R have the following low-T 
expressions: 

I + ^^,1 (^)') - (a + /? (^)' A In a) e-i) f , T « AT,, A « 1 



Fs=i ~ < 



A > 1 



Ri 



+ 0(e-A), A< 1 



2 + C(e-^), A > 1 

where now a, a', (3, jS', and c are constants of order 1, 7j^i = 7t/2 and cjj i = A/tt"^ are the 
coefficients for the specific heat and the susceptibility of the 5 = 1 isotropic model. 

When Tj = and Ta is finite, the system behaves as a single channel Kondo model 
with S = 1, describing at low temperatures a partially screened spin. As Ji is turned on, 
(A ^ 1), the system undergoes another Kondo screening, as can be seen in the temperature 



range Ti <^ T ATj = Ta, and ends up in the infrared as a combination of two screened 
5=1/2 single channel models. The screened S = 1 isotropic system, reached as A = 0, 
differs in essential ways from the above combination. 

We also carried out a numerical solution over the whole range of anisotropy and tem- 
perature (Fig. 1). Note (e.g. in the spin-1/2 case) the two stage quenching of the In 2 high 
temperature entropy when Ta < Ti. For Ta > Ti the quenching occurs in one step at Ta. 
Similarly, the two peaks present in the specific heat for A < 1 merge and move with Ta as 
A is increased past 1. Related observations apply also in the spin-1 case as A is reduced 
past 1. 

We proceed now to interpret our results in terms of fixed-point hamiltonians that capture 
the physics of the model in the infrared limit. Clearly, we reach two different fixed points, 
the isotropic i, and the anisotropic a when we set A = 0, and A = oo, respectively. In the 
screened case the a-fixed point is unstable to turning on Ji {Ti > 0), and flows to the fixed 
point describing two screened spin-1/2 models, to be denoted by g. We conjecture that for 
anisotropy A > 1 the model flows to a line of fixed points connecting g toi, while for smaller 
anisotropy it flows to i. 

In the overscreened case, i is unstable, as remarked earlier and flows to a new fixed 
point g'. We conjecture that for anisotropy A < 1 the model flows to a line of fixed points 
connecting g' to a, while for larger anisotropy it flows to a. In both cases the lines of fixed 
point hamiltonians Hfp[A) describe systems with linear specific heat and zero entropy, but 
not necessarily FL! [^]. 

The structure of the thermodynamic equations reflects the flavor symmetry breaking and 
generalizes to any number of flavors. The various patterns of the S'?7(/)-symmetry breaking 
and their relative strengths will be parametrized by energy scales T„ = Dgn{Ji...Jf), a = 
1, .., /. The scales set the excitation energies and momenta and appear in the thermodynamic 
equations. 
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\nrin = -^e« + G(ln(l + r^^.i) + ln(l + rin+i)), n<f 
\nr]n = G(ln(l + r]n-i) + ln(l + r]n+i)), n> f 

with r^o = 0. As an illustration we discuss the 3-channel problem, with scales Ti, T2, T3. The 
isotropic case is characterized by T3 > 0,Ti = T2 = 0. The low energy physics is NFL for 
5* < 3/2 with Cy ~ T^/^ and zero-temperature entropy S = ln((sin(2S' + l)7r/5)/(sin7r/5). 
When the symmetry is broken to SU{2) x f/(l) a new scale appears. When Ti > 0, T2 = 
the models with 5=1/2 and S = 3/2 show linear T dependence in C^, and constant with 
different coefficients for the Ti > T3 and Ti < T3 regions. The S = 1 case is more interesting, 
since there is a change from linear T dependence to NFL behavior, and a residual T = 
entropy emerges as the parameters of the system are changed from Ti < T3 to T3 < Ti. 
When Ti = 0, T2 > 2 and 5* = 1/2, we have two different NFL fixed points for T2 < T3 and 
T3 < T2 with different values of the residual entropy; the S = 1,3/2 cases have constant 
Cv/T and x- Finally, when all the scales are finite, the system has constant and x fo^^ 
S = 1/2,1,3/2, with the coefficients depending on the relation between the Tj and with 
intermediate temperature regions where the properties of the thermodynamic quantities 
correspond to those described in the previous cases. These considerations generalize for any 
number of flavors /. 

In forthcoming work we study the competition between strength of the couplings and the 
size of symmetry breaking. We also consider the role of the generalized fusion mechanism 
in the 2-impurity Kondo model, and in models (e.g. with 0(A^)— symmetry) giving rise to 
new types of fixed points. Finally, we shall discuss some phenomenological implications. 

We received two preprints UTO) addressing the model by means of Bosonization and the 
Anderson- Yuval Approach. While there is a considerable overlap with our work, there are 
also interesting differences. 

We wish to thank P. Coleman, M. Douglas, A. Gogolin, J. Moreno, and A. Ruckenstein 
for enlightening and stimulating discussions. N.A. wishes to thank P. Woelfle and T. Kopp 
for their warm hospitality at the Institut fuer Theorie der Kondensierten Materie where part 
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of the research was carried out. 



REFERENCES 

[1] See e.g. C.L.Seaman, M.B. Maple, B.W.Lee, S.Ghamaty, M.S. Torikachvili, J.S. Kang, 
L.Z. Liu, J.W. Allen, and D.L. Cox., Phys. Rev. Lett. 67 2882, (1991). 

[2] Ph. Nozieres and A. Blandin, J. Phys. (Paris) 41, 193 (1980). 

[3] N. Andrei and C. Destri, Phys. Rev. Lett. 52, 364 (1984). 

[4] A.M. Tsvehck and P.B. Wiegmann, J. Stat. Phys. 38, 125 (1985). 

[5] P. Schlottman and P.D. Sacramento, Adv. in Phys. 42, 641 (1993). 

[6] L Affleck, A.W.W. Ludwidg, Nucl. Phys. B352, 849 (1991). ibid 360 641, (1991). 
I.Affleck, A.W.W. Ludwig, H.B. Pang and D.L. Cox, Phys. Rev. B 45, 7918 (1991). 

[7] V.J. Emery and S. Kivelson, Phys. Rev. B 47, 10812(1992). A.M. Sengupta and A. 
Georges, Phys. Rev. B 49, 10020 (1994). 

[8] L.A.Takhtajan and L. Faddeev, Russian Math. Surveys 34:5,11 (1979). 

[9] N. Andrei and A. Jerez, Rutgers preprint. 

[10] M. Fabrizio, A. O. Gogolin, and Ph. Noziers, ILL preprints. 



11 



FIGURES 

FIG. 1. (a) Entropy, (b) Specific heat, and (c) Zero field susceptibility for S = 1/2 as a function 

of T/Tj, for different values of the anisotropy parameter A = Ta/Ti. (d), (e) Same as (a), (b), now 
for 5 = 1. (f) Zero field susceptibility at T = as a function of the anisotropy parameter A. 
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